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Abstract 

Post LSND results, sterile neutrinos have drawn attention and motivated the high energy 
physics, astronomy and cosmology to probe physics beyond the standard model consid¬ 
ering minimal 3-1-1 (3 active and 1 sterile) to 3-l-N neutrino schemes. The analytical 
equations for neutrino conversion probabilities are developed in this work for 3-1-1 neu¬ 
trino scheme. Here, we have tried to explore the possible signals of T and CPT violations 
with four flavor neutrino scheme at neutrino factory. Values of sterile parameters consid¬ 
ered in this analysis are taken from two different types of neutrino experiments viz. long 
baseline experiments and reactor-|-atmospheric experiments. In this work golden and 
discovery channels are selected for the investigation of T violation. While observing T 
violation we stipulate that neutrino factory working at 50 GeV energy have the potential 
to observe the T violation signatures for the considered range of baselines(3000 km-7500 
km). The ability of neutrino factory for constraining CPT violation is enhanced with 
increase in energy for normal neutrino mass hierarchy(NH). Neutrino factory with the 
exposure time of 500 kt-yr will be able to capture CPT violation with hcai > 3.6 x 10“^^ 
GeV at 3a level for NH and for IH with hcsi > 4 x 10“^^ GeV at 3a level. 
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1 Introduction 


The standard model of particle physics considers neutrinos to be massless. Sudbury 
Neutrino Observatory n 0 gave evidence of neutrino oscillations which was further 
conhrmed by KamLAND experiment [3]. This landmark research assigned mass to the 
neutrinos and gave a clear indication of new physics beyond the standard model. A sim¬ 
ple stretch in the standard model was able to stand up with the mass of neutrino. In 
neutrino physics the standard three flavour neutrino oscillations can be explained with 
the help of six parameters namely 6 * 12 , 6 * 13 , 6 * 23 , Am^ 2 ) and 5cp- Amongst these six 

parameters, solar parameters( 0 i 2 , Am^ 2 ) atmospheric parameters { 623 , Amh) have 
been measured with high precision. Furthermore, Daya Bay and RENO reactor exper¬ 
iments have strongly constrained the value of mixing angle 613 . Now we are in need of 
such neutrino experiments which can impose tight constraints on the value of 5cp and 
mass hierarchy. Some anomalies popped up while observing appearance channel and dis¬ 
appearance channel of Pe at LSND experiment. While observing z7g appearance 

channel, LSND 000000 was the hrst experiment to publish evidence of a signal 
at Am? ^ leV^. Later in 2002 , MiniBooNE (TU] [H] checked the LSND result for Ug —?• 

(pe — 7 - Ufj) appearance channel. In MiniBooNE experiment, while observing the CCQE 
events rate through Pf.n —>■ e~p{PeP —>■ e+n) above 475 MeV energy, no excess events 
were found but for energies < 475 MeV Pe{^e) excess events were observed. In this way, 
MiniBooNE supported the LSND result. The LEP data [12] [IH| advocates the number 
of weakly interacting light neutrinos, that couple with the Z bosons through electroweak 
interactions, to be 2.984 ± 0.008; thus closing the door for more than three active neutri¬ 
nos. Hence, the heavy neutrino announced by LSND group should be different from these 
three active neutrinos. This higher mass splitting in the standard three active neutrino 
model was accommodated by introducing sterile neutrinos. Sterile neutrinos carry a new 
flavor which can mix up with the other three flavors of standard model but they do not 
couple with W and Z bosons. The number of sterile neutrinos can vary from minimum 
one to any integer N. 

Some cosmological evidences like CMB anisotropies HD US] HE] HZ] m and Big Bang 
nucleosynthesis [12] [20] also stood up with the LSND data. The results reported by the 
combined analysis [21] of Baryonic Acoustic Oscillations(BAO) ‘i7o+PlaSZ-|-Shear-|-RSD’ 
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indicated the presence of sterile neutrinos by stipulating the number of effective neutri¬ 
nos Neff = 3.621q;42, (sterile) = AA^'^^^^eV and giving preference for ANeff = 
Neff — 3.621q; 42 at 1.4 a level and non zero mass of sterile neutrino at 3.4 a level. 
The gallium solar neutrino experiments (gallium anomaly) GALLEX |27j, SAGE [2H! and 
the antineutrino reactor experiments (Pe) hke Bugey-3, Bugey-4, Gosgen, Kransnogark, 
IIL [22] (reactor anomaly) indicated that electron neutrinos and antineutrinos may disap¬ 
pear at short baselines. Such disappearance can be explained by the presence of at least 
one massive neutrino (of the order of 1 eV). Thus, these experiments also indicated the 
presence of sterile neutrino and supported the LSND results. Some constraints imposed 
by the combined fit of reactor, gallium, solar and UgC scattering data are Arrih > leV^ 
and 0.07 < sin^2uee ^ 0.09 at 95% GL [HD|- Eew atmospheric neutrino experiments such 
as IceCube [2I|, MINOS [22| [2S| [Sl|, CCER [33] have also imposed strong constraints 
on sterile parameters. 

The four fiavors of neutrino can be studied in either of the two different neutrino mass 
schemes, 3-1-1 or 2-1-2 schemes [36]. For our work we have selected (3-1-1) four flavor 
neutrino mass scheme. In this framework, Maki-Nakagawa-Sakata (MNS) mixing matrix 
(4x4), includes six mixing angles djj, three dirac phases and three majorana phases. In 
our analysis majorana phases are not taken into consideration. 

Neutrino factory [37] [3H] provides excellent sensitivity to the standard neutrino oscilla¬ 
tion parameters and therefore seems to be one of the promising option to explore and 
reanalyze the global fits for sterile neutrino parameters too. To mention, it provides a 
platform to constrain one of the most searched CP violation in leptonic sector [39] [dfij . 
Hence neutrino factory seems to provide a promising environment for the study of T and 
CPT violation. The neutrino factory set up considered here is based on the International 
Design Study of Neutrino Factory [IDS-NF] [H] [12]. From the measure of APcp we can 
not directly constrain CP phase because the value of APcp in the presence of matter will 
contain in itself some CP odd effects even in the absence of CP phase. Therefore, instead 
of checking APcp, variation in APp can be studied to probe extent of true CP violation. 
We have observed T violation through —>■ Ug golden channel and —)■ discovery 
channel and CPT violation along disappearance channel. 

Our work is organized as follows. In section 2, we illustrate 3-1-1 neutrino matrix parametriza- 
tion. In next section, T violating effects are checked for different channels. In section 4, 
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bounds on CPT violating terms are checked in presence of sterile neutrino. In the last 
section, we have summarized our study and discussed the results observed. 


2 Standard Parametrization in 3 + 1 neutrino scheme 

The 3 (active) + 1 (sterile) neutrino scheme can be looked upon as 3+1 or 2+2 scheme 
depending on the selection of mass ordering of the neutrinos. To check T and CPT viola¬ 
tion we have selected 3+1 scheme for our analysis. In this scheme the flavor eigenstates 
= e,fi,T,s) and mass eigenstates Uj(j = 1,2,3,4) are related by the given unitary 
transformation equation 




(.0 


= u 

i^2 

Vr 


i^3 





( 1 ) 


Here unitary matrix (U) can be parametrized in terms of six mixing angles(di 2 , ^ 13 , 
^ 23 ,^ 14 ,^ 24 ,^ 34 ), three Dirac phases 5i (hi,(52,(^3) and three majorana phases. Majorana 
phases are neglected in our study as they do not affect the neutrino oscillations in any 
realistically observable way. In principle, there are different parametrization schemes for 
the neutrino mixing matrix as their order of sub-rotation is arbitrary. Our selection for 
parametrization of neutrino mixing matrix is 

U = 1 / 34 ( 034 , 0 ) 1 / 24 ( 024 , O)f/i4(0i4, 0 ) 1 / 23 ( 023 , 03 )^/ 13 ( 013 , ^ 2 )^/ 12 ( 012 , 0i) (2) 

A N X N unitary matrix contains N(N-l)/2 mixing angles and (N-l)(N-2)/2 Dirac type CP violating 
phases. It will also contain (N-1) number of additional Majorana Phases if the neutrinos are considered 
as Majorana particles. 
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where Uij{0ij,5i) are the complex rotation matrices in the ij plane, defined as 


[Uij{0ij, Sl)]pg — < 


cosOij 

p = q = i,j 


1 

P = Q7^ hj 



P = i , q = j 

(3) 

-sinOije^^'- 

p = j , q = i 


0 

otherwise 


is arbitrary 

since these matrices commute. When 


neutrinos pass through the earth matter, the charge current interactions (CC) of Ug and 
neutral current interactions (NC) of Ug, Ur with the matter give rise to a CC and NC 
potentials Vg and I 4 respectively. While studying the sterile neutrinos, potential 14 can 
not be neglected. The effective CPT violating hamiltonian (i^/) of neutrinos can be 
expressed as 


^0 

0 

0 

yo 


0 

0 

0 


0 

0 

0 


0 

0 

0 


\ 




^0 

0 

0 


Am|i y 




0 0 0 ^ 

5c2i.2E 0 0 

0 (5c3i.2E 0 

0 0 Sc4i.2Ej 


^AgEAn 0 0 0 ^ 

0 0 0 
0 0 0 

^ 0 0 0 Oy 


(4) 


Here He(n) = 2EVg(n), Vg = 'J^CpNg and I 4 = —GpEn/\/2. Gp is the Fermi constant, 
Ng and A4 are the number density of electrons and neutrons respectively with Ng ~ A4 
in earth matter. The hcp’s are CPT violating terms. Different angular values for unitary 
matrix Uf, are checked in [IHj. In our work we have considered U = Ui,. Hamiltonian Hf 
can be diagonalized to Hp by an unitary matrix U. This can be expressed as 


Hd = U^HfU 


(5) 
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The matrix elements [HD]ii will represent the eigenvalues of Hf. Full analytical ex¬ 
pressions for neutrino oscillation probabilities are developed in this work by using time 
independent perturbation theory. In an attempt to apply perturbation we have defined 
few oscillation parameters in terms of perturbative parameter 77 , where 77 = 0.18 . The 
neutrino oscillation parameters can be rewritten as 
Ou = XuV 
^24 = X 24 ?? 

^34 = XUV 

0l3 = Xl3V 

^23 — ^23 ~ 1 / \/2 = X 23 V 


We treat 


Amii 


Oijf). Now the Hamiltonian Hf can be written as 


Hf = 


2E 


[Ho + H^rJ + H2v" + Oiv^)] 


( 6 ) 


where Hq, Hi and H 2 are the hamiltonians corresponding to zeroth, first and second order 
in Tj respectively. The evolution equation for neutrino oscillation probability is defined as 


Pa/S =\ SiSa{t,to) \'^ (7) 

where S{t,to) is the evolution matrix of neutrino which is also called oscillation proba¬ 
bility amplitude 

I i/(t) >= I > ( 8 ) 

The evolution matrix of neutrinos in terms of eigenvalues of H^ can be written as 

4 

S^a{t,to) = (9) 

i=l 

where L = t — to 

From equation (7) the neutrino oscillation probability Pa /3 from flavor a to flavor /3 can 
be written as 

4 

Pa^ =1 P (10) 

i=l 

This is the general form of equation for neutrino oscillation probability. 
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3 T violation in (3+1) framework 


In neutrino oscillations the flavor conversion probabilities from flavor a to flavor (5 can 
be written as 


Pva—yyp — ^ij 4 ^ ^ Re 




sin^ Ajj + 2 Im 


Ua^K,u;,u^, 


sin2Ajj (11) 


i>j 


i>j 


Redefining the above probability equation as sum of Pcp-even and Pcp-odd terms 


-Pa/3 P{va^vi}) PcP—eveni^a t Z^^) + PcP—oddiPa t 


( 12 ) 


CP even terms are CP conserving and can be written as 


~ ^ij 4 ^ ^ Re(JJ^-Uj^j) sin A^j 

i>j 


(13) 


CP odd term are CP violating and can be written as 


P{va^up) — P{i/a^ui 3 ) — 2 ^ ) Im{UaiU^fi^jUi 3 j) sin 2Ajj 

i>j 


(14) 


Assuming CPT to be conserved, the magnitude of CP violation (APcp) will be equal to 
the magnitude of T violation (APp), i.e. 


APcp I—I AR 


T 


(15) 


therefore we can write 






(16) 


When neutrinos passes through the earth matter, the interaction of neutrinos with matter 
gives rise to an extra potential. This potential is positive for neutrinos and negative for 
antineutrinos leading to different eigenvalues of hamiltonian for them. Further, this 
difference in hamiltonian for i^’s and P’s give rise to fake (extrinsic) CP violation. Hence, 
check on T violation appears to be a better choice in the presence of matter. From 
equation (16) the T violation can be looked upon as 
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( 17 ) 


(APT)a /3 = ^ = 4 Im{UaiUlfilfijjj) sin 2Ajj 

i>j 

If we consider llajU^j = and Ajj = 2AEjk/L = 2{Ej — Ek)/L the above equation 
becomes 


(APr)a /3 = 4 /m(I//V/“’) sin (AP.^L) (18) 

j<k 

~ 13a ~ /3a* 

The term Im{Vj 14 ) is known as Jarlskog factor and EjS are energy eigenvalues of 

hamiltonian in matter. 

{APt)^^ = 4 /m(I//“l/ 2 ^“*) sin {AE^L) + 4 Jm(I//“l/ 3 ^“*) sin (AP 13 L) 

+ 4 /m(I/ 2 ^“l 4 ^“*) sin {AE23L) + 4Im{v/'^v/'^*) sin (AEuL) (19) 
+ 4 /m(I/ 2 ^“l 4 ^“*) sin (AP 24 ^) + sin (AP 34 L) 

The energy eigenvalues in matter can be connected to the energy eigenvalues in vacuum 

~ ^ ~ ^ ~ a/3 ~ a/3 

by the relations Pi = AP 31 , E 2 = 0, P 3 = Ag and P 4 = AP 41 . The terms I 4 ,14 , 

14 °^^, 14 °^^ in matter can be calculated with the help of the terms 14 ^^, in 


vacuum (where = UajU^j) using the following expressions 

14 ^“ = -AEA^AEA^{Vf^E2E3 + {E 2 + E3)R^^ + (20) 

14^“ = AEA^AE^^^iVtEsEi + (P3 + (21) 

14 ^“ = -AP3-/AP3441//“PiP2 + (^1 + P2)P^“ + (22) 

1?/“ = l//“ (23) 

where 

= {Aivr + 14^^/2) - A^a - + APsiVf“ + AP2iP2^“ (24) 
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= vt{Al^ + - A(A^^ + A^^)(vr + Vfy2)} 

- AE3i(AB3i + + A^^)Vf^ - AB2i(AB2i + + A^^)Vt 

( 25 ) 

+ AABaiiVtVi^ + VtVi^ + + V^Vl^) 

+ AAE2i{vtvi^ + + vtvi^ + if *V7“) 

Aa = AeSae—AnSas IS the diagonal element of the matter potential matrix in four neutrino 
scheme and A is the diagonal element of matter potential matrix in three neutrino scheme. 
Since T violating effects can only be studied in appearance channels so a 7 ^ /3. In an 
effort to put constraints on APt we have studied two appearance channels. These are 
t (golden channel) and —>■ (discovery channel). 

The T violation probability difference expression for the golden channel can be expressed 
as 


(APr)Me = 


-4 


{Amy2E)^^{^-Aml/2E) 


X 


^ S^%s^,c^,S 24 ^AmlJ 2 E + (^ + AmlJ 2 E)R^^ + smAml,L/ 2 E 


+ 


[-AmlJ2E)^(^ - AmlJ2Ey- 


X 


[ A ^ S‘>‘][(AmlJ2E)W‘^ + 5'"] sin (Am§,/2B - T)l+ 


-A™y2£)(i 1 A™i./2£)(i)^ [suCuS^l^^mlJ2E + + A™L/2S)iJ- 


^2 


A. 


+ X [{Aml-^/ 2 E)R^>^ + sin -^L 

2 E 


( 26 ) 


Since large value of Amh gives rise to rapid oscillations, hence Amh terms can be 
averaged out. Solving the above expression up to the power sfj we get 


{APT)^e = 4C13C14C24S13S23S14S24 sin {82 - S3) 


A. 


(Ae — A 31 ) 


sin A'm\iL/2E 


+ 4C13C44C24S13S23S14S24 sin (^2 - ^3) 


A 31 


( 27 ) 


^ ^ ^ - sin 2 Ae 

Ae(Ae - A 31 ) 
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( 28 ) 


{APT)i^e = 4 C 13 C 14 C 24 S 13 S 23 S 14 S 24 sin ((^2 - (^ 3 )[ 


(Ae — A 31 ) 


sin Am\iL/2E 


+ 


A 2 

^31 


Ae(Ae Asi^ 


sin 2 Ap 


Ag = A^LjAE is matter dependent term. The change in Ag will change the value of 
(APT)Me- 

Further we have developed equation of AP^ for discovery channel. The discovery chan¬ 
nel is not very useful in the standard three neutrino flavor framework, nevertheless while 
studying physics beyond three active neutrino flavor framework, it becomes very impor¬ 
tant. For discovery channel(z 2 ^ Vr) the probability difference is given as 

-4 


(APt),., = 


(Aml/2E)^M^ - AmlJ2E) 


X 


2E\2E 
2 


+ <S'^''][cMC24S24S34^Am3i/2P Aml^/2E)R'^^ -F sin Am 3 iL/ 2 P 


+ 


AmlJ2E)^{^-AmlJ2Ey 


■X 


+ S-f][{AmlJ2E)R'f + V'-] sin (AmlJ2E - A)l+ 
(-An,y2P)(^-An.l./2P)(^)4 -tc.s.S34An.L/2P P (A + AmlJ2E)R^^ 


^2E 

2 


Ag 


+ X [(Am^ 4 / 2 P)P^^ + sin 


Solving the above expression up to the power sh we get, 


(29) 


{APT)^r = 4Ci3cf4C23C24C34S23S24S34 sfo siu Aml-^L/2E 

^31 


4 C 43 C 14 C 23 C 24 C 34 S 23 S 24 S 34 sin (^ 3 )-^ sin 2 A, 


31 


(30) 


Ag 


iAPT)f,r = 4cf3C44C23C24C34S23S24S34 sfo ((^3)[x^ sin Aml-^L/2E 

^31 

A 3 I . o A 1 

- — sm 2 Ag] 


( 31 ) 


The /\Pt for three neutrino framework is given by 
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A(Pr)3x3 4C12C?3C23S12S13S23 sin (S 


(32) 


Keeping the best fit values of neutrino oscillation parameters and assigning maximum 
value to dirac phases i.e keeping mod of sin of dirac phases to be unity will lead us to 
maximum value of APt- This assumption will render the maximum limit on the bounds 
which can be imposed on T violation arising due to the presence of dirac phases if all 
other oscillation parameters are known with utmost accuracy. From equation (32) gives 
the value of (APr)rnaa; = 0.137 for three neutrino flavor framework [IH]- This value is 
independent of the selection of probing channel and presence of matter effects. Whereas 
in 4 flavor framework it will depend on the selection of channel through which we want 
to probe CP or T violation and it will vary with matter effects too. Within 4 flavor 
neutrino framework the magnitude of APt will depend on active flavor neutrino mixing 
angles (known with accuracy), sterile neutrino mixing angles (still needs better bounds), 
matter effects, baseline, energy and dirac phases (not known) . Imposition of constraints 
on dirac phase (three neutrino flavor) or phases (four neutrino flavor) is still in research 
phase. 


0.5 


Q. 

< 

-0.5 

-1 

-1.5 
3i 

(£ 

Figure 1: The variation in APx with baseline for two different appearance channels 

The probability differences {APt) 4 x 4 in neutrino sector are represented by equations 
(28) and (31) for 12 ^ —>■ Ug and 12 ^ Ur channels respectively. The values of sterile 
parameters used in the above mentioned equations are taken from 


Long baseline constraint, 4 GeV 
Long baseline constraint, 50 GeV 
Reactor constraint, 4GeV - 
Reactor constraint, 50 GeV - 
3 flavor 



100 3500 4000 4500 5000 5500 6000 6500 7000 7500 

L(Km) 



,)Variation in (APr)//e along the baseline (b)Variation in {APT)^ir along the baseline 
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(i) long baseline neutrino oscillation experiments jlTj 
0U < 6.7°, 024 < 3.3° and 034 < 6.3° 

(ii) short baseline reactor and atmospheric experiments [IHj [IH] 

014 < 10°, 024 < 12° and 034 < 28° 

The variations in (APr)^e and {APt)^t are checked along the baseline for two different 
energies i.e 4 GeV and 50 GeV. Plot (a) of Figure 1 reflects a very small variation in 
(APr) 4 x 4 in comparison to (APr) 3 x 3 for golden channel when it is checked with two 
different energies and two sets of sterile parameter values. Plot (b) of Figure 1 reflects a 
reasonable variation in (APt) 4 x 4 in comparison to {APt) 3 x 3 for discovery channel when it 
is checked for 50 GeV energy and sterile parameters are taken from reactor+atmospheric 
experiments. The Pr comes up as the promising channel to observe the signatures 

of T violation. From the analysis we conclude that neutrino factory operating at 50 GeV 
has the potential to capture the signatures of T violation through p^ —Pr channel if true 
sterile parameter values are equal to that taken from reactor+atmospheric experiments. 
At the same time if the upcoming neutrino experimental setups captures APt value 
above 0.137 then we can stipulate presence of some new physics beyond three active 
flavor neutrino physics which is responsible for the enlargement observed in APt value. 


4 CPT violation in 3 + 1 scheme 

GPT invariance is one of the most fundamental symmetries of nature. GPT conserva¬ 
tion indicates the invariance in the properties of physical quantities under the discrete 
transformations such as charge conjugation (G), parity inversion (P) and time reversal 
(T) along with the invariance under lorentz transformation. GPT invariance is one of 
the symmetries of local quantum field theory which implies that there is an important 
relation between GPT invariance and Lorentz invariance. If GPT invariance is violated, 
Lorentz invariance must violate but if Lorentz invariance is violated it is not necessary 
that GPT invariance must violate. In our work Pj^ —>• Pj^ disappearance channel is probed 
to check GPT violation. The GPT violating probability difference can be written as 

A PSr = P„0 - Pfs (33) 
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The intrinsic CPT violation arises due to the violation of CPT invariance theorem. A 
hamiltonian Hf containing CPT violating terms is defined by equation ( 4 ) and the gen¬ 
eral form of neutrino oscillation probability is mentioned in equation ( 10 ). The terms 
Uai,Uaj^^li ^l3j expression ( 10 ) are the elements of the unitary matrix U. 

The construction of unitary matrix U with the help of eigenvalues and eigenvectors of 
hamiltonian Ho,Hi and H2 is mentioned in the Appendix. After the formation of unitary 
matrix U, we have developed the neutrino oscillation probability equations up to second 
order in rj. Since Amh is large, so we average out the effects produced due to Am^i in the 
probability equations. Neutrino oscillation probabilities for different oscillation channels 
containing CPT violating parameters can be developed as 


Pee = l 


2 ^ 14 ' 


4di3^(A3i + hc 3 iL/ 2 )^ 


sin^ (A 31 -I- (5c3iL/ 2 - Ae) 

(A 3 l+( 5 c 3 lL/ 2 -Ae )2 


( 34 ) 


P 

e/j, 


= P = 

er 


2 ^ 13(^31 + SC 31 L/ 2 Y 


sin^ (A31 -F 5csiL/2 - Ae) 
(A31 + 5 c^iL /2 — Ae)^ 


( 35 ) 


= 1 - 26'24 cos^ (A31 SC31L/2) - (1 - 89 ^^)sin‘^{A31 SC31L/2) 

A (^12^12 ~ 2^24^34 cos (53A„) sin2(A3i + (5 c3iL/2) 

+ / [^(^31 + SC31L/2) sin Ae cos {A31 + SC31L/2) 

(A31 + ()C3iL/2 - AeY 

sin (A31 + 5 c^iL /2 — Ag) — (A31 + 5 c^iL /2 — Ae)AeSin2A3i + ScsiL/ 2 ] 


( 36 ) 


PfiT — sin^ (A31 -|- dc3iL/2) — { 89 ^^ + 024 + ^34) (^ 3 i + SC31L/2) 

~ (Ci2^12 + 2024034 COS S^An) siu 2(A3 i + 8c^iL/2) — - fj^ - ^/\ \2 

(A 31 + ()C 3 iL /2 - AeY 

[ 2 (A 3 i + 5c^\L/2) sin (A 31 + 5cz\L/2) cos Ag sin (A 31 + 5c^\L/2 — Ag) 

- (A 31 SC 31 L /2 - Ae)AeSin 2 (A 3 i 5c3iL/2)] 

+ 024034 sin 03 sin 2(A3 i + 0 C3iiv/2) 


Pf,s = 2024 + (034 - 024 ) sin^( 2 ^ 3 i + 0C3iL/2) - 024034 siu 03 siu 2(A3i 0C3iL/2) ( 38 ) 
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For small angles {9ij ~ sin6*jj ~ Sij) these oscillation probabilities can be written as 


1 - 2514^ 


~ 4si3^(A3i + 5c3iL/2)^ 


sin^ (A 31 + (5c3iL/ 2 - Ae) 
(A3l+hc3lL/2- Ae)2 


(39) 


p = p 

^ e/j, ^ er 


2'5i3(^31 + ^C 3 iL/ 2 )^ 


sin^ (A 31 + 5c^iL/2 — Ag) 
(A 31 + 5c^iL/2 — Ag)^ 


(40) 


PiJLiJL — 1 ~ 2^24 cos^ (A31 + bc‘^\Lj 2 }j — (1 — 85^3) sin^ (A31 + bc‘^\Lj 2 }j 

+ (C 42 A 12 - 2524534 cos 53A„) sin 2 (A 3 i + 5 c3iL/2) + /\2 

(A31 + dCziL /2 - I\e) 

X [ 2 (A 3 i + SC31L/2) sin Ae cos (A 31 + 5c^\Ll2^ sin (A 31 + 5c^\Ll2 — Ag) — 
(A 31 + ( 5 c 3 iL /2 - Ae)Ae sin 2 (A 3 i + 5 c 3 iL/ 2 )] 


(41) 


P^r = sin^ (A31 + hc3iL/2) - (85^3 + 5^4 + 5I4) sin^ (A31 + 5 c^iL/ 2 )- 

(C42 Ai2 + 2524534 COS 53A„) sin 2(A3 i + 5c3iL/2) - /\2 

(A31 + dC3iL/2 - Ae)^ 

X [2(A3i + 5c3iL/ 2) sin (A31 + 5c3iL/2) cos Ae sin (A31 + 5 cziL /2 — Ag)- 

(A31 + 5 c^\L /2 — Ag) Ag sin 2(A3i + 5 c^iL/ 2 )\ + S24534 sin ^3 sin 2(A3i + 5 c^\L/ 2 ) 

( 42 ) 


= 2524 + ('sL - ■sL) (A 31 + ^C 3 iL/ 2 ) - 524534 sin ^3 sin 2(A3i + 5c3iL/2) (43) 

In order to analyse CPT violation at probability level in 4 flavor neutrino framework the 
value of considered in our analysis is given by 

A PaP = [{Pal3)‘iA5cij^0 “ [(-PQ/3)4i.]5cij=0 (44) 
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(c) 


(d) 


Figure 2: The oscillographs (a) and (b) demonstrates the variation of with base¬ 

line and energy for normal hierarchy whereas (c) and (d) illustrates the same for inverted 
hierarchy. For the oscillographs (a) and (c) the values of sterile parameters are selected 
from long baseline experiments while for (b) and (d) the values of sterile parameters are 
taken from reactor and atmospheric experiments. 
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(a) Normal hierarchy (b) Inverted hierarchy 


Figure 3: The oscillographs demonstrate the variation of with baseline and energy. 

The values of sterile parameters are selected from reactor+atm experiments. Left and 
right oscillographs are for normal and inverted hierarchies respectively. 

Neutrino factory setup considered for analysing CPT violation is taken from the refer¬ 
ences [HO] [SI] [S2| jSS] and jSl]. The experimental setup and detector specifications con¬ 
sidered in our analysis are mentioned below. Neutrino factory setup consist of 1.4 x 10^^ 
useful muon decays per polarity, with parent muon energy = bOGeV . We have 
done our analysis for 10 years running of neutrino factory. In particle physics meaningful 
observations always demands a detector with very good energy and angular resolutions. 
This view point lead us to select Liquid Argon detector for particle detection. The energy 
resolution of the detector for muon is a{GeV) = 0.20/■ 

A near detector is placed at a distance 20 m from the end of the decay straight of the 
muon storage ring. Effective baseline (Le//) is used in place of baseline (L), which is 
calculated using L^ff = d{d s) [S3]. Fiducial mass of near detector is 200 tons. Pres¬ 
ence of near detector will minimize the systematic uncertainties in our observations. 

A 50 Kt far detector is placed at a distance of 7500 Km. The systematic uncertainties 
considered for this analysis are given in Table 1. An uncertainty of 5% on matter den¬ 
sity [SH] [SI] is also considered in our work . The simulated environment of the neutrino 
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factory is created with the help of GLoBES [58j [59j. The analytical equations for four 
flavor nentrino conversion probabilities derived in this work are defined in the probability 
engine of the software. The best fit values of oscillation parameters m ra are men¬ 
tioned in Table 2. Sterile parameter values mentioned in the Table 2 represents best fit 
values for =0.1 eE^. 

Table 1: Systematics 


Systematic uncertainties values 

Flux normalization 2% 

Fiducial mass errors for near detector 0.6 % 
Fiducial mass errors for far detector 0.6 % 
energy calibration error for near detector 0.5 % 
energy calibration error for far detector 0.5 

shape error 10 % 

Backgrounds 10“^ 


Table 2: Best fit values of the oscillation parameters 


Parameter 

Best fit values 

^12 

34.4° 

^13 

8.50° 

6*23 

45.0° 

6*14 

o 

O 

o 

^24 

3.3° ; 12° 

6*34 

6.3° ; 28° 

Am2i^ 

8 X lO-^eF^ 

Amsi^ 

2.5 X 10“^eF^ 


The constraints on CPT violating parameters 5 c2i and within two and three neu¬ 
trino frameworks are mentioned in references |46| [F)T] [H2] jHH] |64| |65| and [HS! • In present 
work we are trying to check the neutrino factory potential to capture CPT violating sig¬ 
natures in presence of sterile neutrino. As the mass hierarchy determination is yet in 
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research phase therefore in an attempt to make this work relevant we have analysed CPT 
invariance for both the mass hierarchies. Initially CPT violating signatures are checked 
at probability level. The value of is estimated by substituting equation (41) in 

equation (44) for the channel Vn ^ Vn- The variation in AP^^^ and AP^P^ with baseline 
and energy is shown in Figure 2 and Figure 3 oscillographs respectively. The total CPT 
violation captured by any experiment will be the sum of extrinsic CPT violation (CPT 
violation arising due to matter effects) and intrinsic or genuine CPT violation (which we 
are probing in present work). In an endeavour to constraint intrinsic CPT violating 
parameters we must look for places where extrinsic CPT violation is negligible or very 
less. With three active neutrinos the extrinsic CPT violation is checked in reference Ea 
whereas with 3 (active) + 1 (sterile) neutrinos it is checked in reference [HE]- These 
references conclude that extrinsic CPT violation for energies 4 GeV- 6 GeV is negligible 
at shorter baselines, roughly less than 2000 km. Whereas for longer baselines this effect 
decreases with energy. Equation (44) of our work will check the presence of pure CPT 
violation arising in the presence of sterile neutrino at probability level. The values of CPT 
violating parameters considered while plotting oscillographs are Sc^i = 4 x 10“^^ GeV 
and (5c2i = 3 x 10“^^ GeV. Looking at normal and inverted hierarchy oscillographs (Figure 
2) we can observe the presence of pure CPT violating signatures at shorter baselines i.e. 
from 1300 km-2000 km for 4 GeV-6 GeV energies. The references [SZ! [SB], which speak 
about extrinsic GPT violation have recorded very weak or almost negligible signatures of 
extrinsic CPT violation at the above mentioned energies and baselines. Hence baselines 
from 1300 km-2000 km with neutrino energies in the range 4 GeV-6 GeV are favourable 
for probing GPT violation with neutrino factory. In Figure 2 while looking at normal 
hierarchy oscillographs we observe AP^J’^ =-0.05 along baselines 4000 km-7500 km for 
energies 12 GeV-30 GeV. The inverted hierarchy oscillographs of the same Figure cap¬ 
tures AP^^^=0.25 and 1 for sterile parameters taken from long baseline experiments 
and reactor-|-atmospheric experiments respectively. This probability difference can be 
observed for baselines 4000 km-7500 km and energies 20 GeV-50 GeV. 

After examining the presence of pure GPT violation at probability level we go ahead to 
observe the signatures of the same with realistic proposed neutrino experiments i.e. neu¬ 
trino factory. The specifications of neutrino factory considered in our work are mentioned 
earlier. Liquid argon detector seems a reasonable choice to grab the signatures of leptons 
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in the considered energy range. The rate(event) level analysis depends on mathematical 
formulation( oscillation probability), physics(types of interactions) and R & D (source 
properties and detector properties) of the experiment. Looking at equations (34) to (44) 
we found that hc2i appears with Am^i term and Sc^i term appears with Amh term. The 
solar and atmospheric mass square difference {Am^i and Amh) are of the order of 10“^^ 
and 10“^^ respectively. Hence, any change in mass terms due to the presence of CRT 
violating parameter will be better observed in term. 

To hook CPT violating impression with neutrino factory we have investigated some ob¬ 
servable parameters like R, AR and asymmetry factor. These terms are defined by 
equations (45) and (46). The ratio R and ratio difference AR are examined as 

^ (46) 

Niyij, 

where u^) denotes number of muon neutrinos reaching at detector as muon 

neutrinos and producing a lepton and A^(i^ —>■ z^) denotes number of anti muon 
neutrinos reaching at detector as anti muon neutrinos and producing a lepton. In 
AR, {R 4 ,u)scij^o denotes the ratio R in presence of CPT violating terms and {R 4 ,u)scij=o 
denotes the ratio R in absence of CPT violating terms. 
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Figure 4: The variation of R with energy in the energy range 4-15 GeV. Values of sterile 
parameters considered in plot (a) and plot (b) are selected from long baseline experiments 
and reactor-|-atm experiments respectively. These observations are made for different 
values of CRT violating parameter (i) 5c^i = 0 (setting OPT violating parameter to 
zero) (ii) hcai = 3 x 10“^^ GeV (hi) = 3.5 x 10“^^ GeV (iv) hcai = 4 x 10“^^ GeV. 
For all the observations 8 c 2 i = 3 x 10“^^ GeV 




(a) (b) 

Figure 5: The variation of AR with energy in the energy range 4-15 GeV. Values of sterile 
parameters considered in plot (a) and plot (b) are selected from long baseline experiments 
and reactor-hatm experiments respectively. These observations are made for different 
values of CPT violating parameter Sc^i; (i) Sc^i = 0 (setting CPT violating parameter to 
zero) (ii) hc3i = 3 x 10“^^ GeV (hi) 5c^i = 3.5 x 10“^^ GeV (iv) 5c^i = 4 x 10“^^ GeV. 
For all the observations hc2i = 3 x 10“^^ GeV 
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5c3i=0,NH - 

5c3i=3x10"GeV,NH 

Scgn =3.5x10 "GeV,NH - 

5c3i=4x10’^^GeV,NH- 

5c3i=0,IH 

5c3H=3x10“GeV,IH 

5c3i =3.5x1 O^GeVJH - 

5c3i=4x10'^^GeV,IH- 




(a) (b) 


Figure 6: The variation of R with energy in the energy range 15-50 GeV. Valnes of sterile 
parameters considered in plot (a) and plot(b) are selected from long baseline experiments 
and reactor+atm experiments respectively. These observations are made for different 
values of CRT violating parameter Sc^i (i) Scsi = 0 (setting CPT violating parameter to 
zero) (ii) (5c3i = 3 x 10“^^ GeV (hi) Sc^i = 3.5 x 10“^^ GeV (iv) Sc^i = 4 x 10“^^ GeV. 
For all the observations Sc 2 i = 3 x 10“^^ GeV 


In presence of matter the observable R will not be equal to one, even if pnre GPT 
violation is absent. It will be equal to a numerical value representing the ratio of neutrino 
and antineutrino interaction cross-sections. If we want to analyse the extent of deviation 
produced by pure GPT violation, we have to hide or filter ont the deviation prodnced by 
any other phenomenon. In an attempt to filter out pure CPT violating contribution from 
the total observed deviation we take into record a new observable AR. This parameter is 
defined in the eqnation (45). Figures 4 and 6 demonstrate the variation in R with energy 
whereas Fignres 5 and 7 exhibit variation in AR with energy for baseline 7500 km. We 
observe that at long baselines pure CPT violating effects get smaller with increase in 
energy. The presence of CPT violation signatures can be observed with neutrino factory 
and it can be checked by looking R and AR plots (Figure 4- Figure 7) for different valnes 
of CPT violating parameter Sc^i. As we know that in presence of sterile neutrino the 
manifestation of pure CPT signatures depends on the values of sterile parameters, hence 
the entire analysis is performed with two sets of best fit valnes of sterile parameters which 
were examined by different nentrino experiments. The results from neutrino factory with 
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(a) (b) 

Figure 7: The variation of AR with energy in the energy range 15-50 GeV. Values of sterile 
parameters considered in plot (a) and plot (b) are selected from long baseline experiments 
and reactor+atm experiments respectively . These observations are made for different 
values of CPT violating parameter Sc^i (i) Scsi = 0 (setting CPT violating parameter to 
zero) (ii) (5c3i = 3 x 10“^^ GeV (iff) Sc^i = 3.5 x 10“^^ GeV (iv) Sc^i = 4 x 10“^^ GeV. 
For all the observations Sc 2 i = 3 x 10“^^ GeV 


sterile parameter values obtained from reactor+atmospheric experiments exhibit larger 
deviation in observables R and AR in comparison to the results obtained with sterile 
parameter valnes taken from long baseline experiments. These observables are checked 
for both mass hierarchies. From the Fignres 4,5,6 and 7 we comprehend that after 15 
Gev there is a flip in sign of the observables for both the hierarchies. At the same time 
the amount of deviation measnred for pure GPT violating effects are different for NH and 
IH for the same energy and baseline. 

The next observable asymmetry factor is defined as 

A (E) = ^ ^ i^;.y"(E) 

' N(iy„ -> !'„)”“'■(£;) -> 

This ratio is determined by using far and near detectors. The variations in asymmetry 
factor with energy for baseline 7500 km are shown by Figures 8 and 9. These figures reflect 
the variations in observable for different values of Sc^i (CPT violating parameter) and 
for both mass hierarchies. The hc3i=0 will reflect values without any contribution 
from CPT violating terms. The asymmetry factor increases with the increase in the 
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value of Scsi. An enhancement in magnitude of asymmetry factor is also observed with 
the increase in values of sterile angles. Hence more stringent bounds on sterile parameters 
are required to check the extent of CPT violation. 




(a) (b) 

Figure 8: The variation in asymmetry factor A^(E) as a function of energy in energy 
range 4-15 GeV. Values of sterile parameters considered in plot(a) and plot(b) are taken 
from long baseline experiments and reactor-|-atm experiments respectively 
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E(GeV) 



(a) 


(b) 


Figure 9: The variation in asymmetry factor A^(E) as a function of energy in energy 
range 15-50 GeV. Values of sterile parameters considered in plot(a) and plot(b) are taken 
from long baseline experiments and reactor-|-atm experiments respectively 
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8c3i (GeV) 




8C31 (GeV) 

Figure 10: The contours are plotted in hcai -Sc 2 i plane for 90% C.L. with CPT vi¬ 
olating terms Scsi = Sc 2 i = 3.0 x 10“^^ GeV taking energies 15 GeV,25 GeV and 
50 GeV respectively. In top two plots sterile parameters values are taken from reac- 
tor-|-atmospheric experiments while for the bottom plot these values are taken from both 
the reactor-|-atmospheric and long baseline^^periments. 






8c3i (GeV) 




8C31 (GeV) 

Figure 11: The contours are plotted in ^031 -Sc 2 i plane for 90% C.L. with CPT violating 
terms Scsi = 4.0 x 10 “^^ and Sc 2 i = 3.0 x 10 “^^ GeV taking energies 15 GeV,25 GeV 
and 50 GeV respectively. In top two plots sterile parameters values are taken from 
reactor+atmospheric experiments while for the bottom plot these values are taken from 
both the reactor+atmospheric and long ba^Gine experiments. 
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Figure 12: as a function of Sc^i is shown. The green curve will represents inverted 

hierarchy and red curve represents normal hierarchy as true hierarchy. 

In our work we have imposed bounds on CPT violating parameters 5c^i and 5c2i at 
90% C.L. Figures 10 and 11 demonstrate contours in 5c^i and 5c2i plane with true value 
of CPT violating parameters as 5c3i = 3 x 10“^^ GeV,5c2i = 3 x 10“^^ GeV and (^cai = 
4 X 10“^^ GeV,(^C2i = 3 X 10“^^ GeV respectively. Each Figure consists of three plots at 
three different energies 15 GeV, 25 GeV and 50 GeV for baseline 7500 km. These plots 
illustrate bounds on CPT violating parameters at mentioned energies. The selection 
of three different energies are based on the results of previous observations (i.e R, Ai? 
and A^). At selected energy 15 GeV, change in sign (+ve to -ve) is observed in the 
observables while studying effects of CPT violation for considered baselines and energies. 
This observation makes 15 GeV energy important for studying CPT violating effects. 
A proposal of neutrino factory producing neutrino beam of 25 GeV muons is described 
in reference j69] whereas in a different proposal we have a 50 GeV muon beam for the 
production of neutrinos in neutrino factory ng. Therefore, by checking at the extent 
of bounds imposed on CPT violating parameters with energies 25 GeV and 50 GeV we 
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want to check that by what order the results will improve if we move towards higher 
energies. By looking at different energy contours we conclude that amongst the three 
selected energies, 50 GeV energy is the best suited energy to constrain OPT violating 
parameter hcai, if nature allows NH to be true hierarchy. At the same time we observe 
that for long baseline experiment IH will be favourable hierarchy for determination of 
bounds on CPT violating parameters for energies less than 15 GeV. 

As discussed earlier that, out of two parameters 5c^i and 5c2i considered in our analysis, 
the variation in Scsi will produce larger variation in the detectable observables which are 
used in our work for checking CPT violation. Figure 12 shows value of as a function of 
CPT violating parameter It is plotted by marginalizing over oscillation parameters 
in Sa range of their best fit values and 5cp from 0 to 27r. Looking at figure we 
observe that the presence of CPT violation can be detected for ^ 3.6 x 10“^^ GeV 
with neutrino factory for NH within 3cr limit. 


5 Conclusions 

Neutrino factory will provide us a potential setup for observing T violation and setting 
significant bounds on CPT violation in neutrino sector. In four(3+l) neutrino flavor 
framework the angular mixing parameters of three active neutrinos are well constrained 
while the sterile parameters still needs better bounds on them. With the change in the 
value of sterile parameters a notable variation in bounds on CPT violating parameter 
and on the extent of T violation is captured by neutrino factory. Hence, well constrained 
values of sterile parameters will allow any neutrino experiment to impose better con¬ 
straints on T violation and CPT violating parameters. Amongst two selected sets of 
values of sterile parameters i.e. from long baseline experiments and reactor-|-atmospheric 
experiments we observed that neutrino factory potential for investigating T and CPT 
violation enhances when the sterile parameters values are equal to those which are con¬ 
strained by reactor and atmospheric experiments. Neutrino factory with 50 GeV energy 
is sensitive to probe T violation when true values of sterile parameters will be equal to 
those predicted by reactor-|-atmospheric experiments. We stipulate that a pure CPT 
violating effects can be observed along short baseline i.e 1300 km-2000 km with energies 
4 GeV to 6 GeV where extrinsic CPT violation is negligible. On the other hand at long 
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baselines we can observe these effects with energies in the range 20 GeV- 40 GeV along 
baselines 4000 km-7500 km. CPT violating parameters ^ 3.6 x 10“^^ GeV for NH 
and 5c^i ^ 4 x 10“^^ GeV for IH will make neutrino factory capable to capture signatures 
of CPT violation at 3 a level. 


A Eigenvalues and Eigenvectors of Hamiltonian to 
second order 

Using the time independent perturbation theory we calculate the eigenvalues and eigen¬ 
vectors of hamiltonian Hf up to the order of correctly. 

Eigenvalues of Hq are given by 


= tte + ttn, = an, E^^’ = Gg + 1, El""' = a 


( 0 ) 


ho) 


where a = 


Aml^ 


and a^^n 


A 

^e,n 


Eigenvectors of Hq are given by 


(A.l) 


1 


0 


0 


0 

0 

II 

—C23 

II 

-S 23 

0 

II 

0 

0 


S 23 


C 23 


0 

0 


0 


0 


1 


(A.2) 


Eigenvalues and eigenvectors for Hi are calculated by using equations (A.3) and (A.4) 
respectively. 


'(1) =< U.0 1 

J J 1 

H,\Vf> 

(A.3) 

^ E1 > 

<Vf\H,\V°> 

(A.4) 

eP - 
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Eigenvalues and eigenvectors of H 2 can be calculated with the help of zeroth and hrst 
order eigenvalues and eigenvectors mentioned in equations (A.5) and (A.6). given as 


=< V«\H,\V}> 


(A.5) 


>= 

3 


3 


where 






k^j 


vr > 


E 

fc'/j 


n I I n ht n 


^kj^k'j 


kj ^jj 


Vkj =< «■“ I r I / > and Eij = Et-lfl 

The total eigenvalues and eigenvectors of up to second order is given by 


(A.6) 


Etotai = + E] + E? 


3 ' 3 ' 3 


(A.7) 


Vtotai = V3 + V3 + V3 


(A.8) 


Using the set of four normalized eigenvectors we form the unitary matrix U as. 


U = 



(Vl^h 

{Vlm)3 

(Ui^)4 

{V2m)l 

iV 2 m )2 

{V2m)3 

(U2^)4 

{V3m)l 

(V,n)2 

{V3m)3 

iV3m)4 

{V 4 J 1 

iV4m)2 

{V 4 J 3 

(h4m)4 


(A.9) 


where Vjm is normalized vector. 


Now hamiltonian Hf can be diagonalised by using the above derived unitary matrix U 
and the diagonalized hamiltonian Hd can be expressed as 

Hd = U^HfU (A.IO) 
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